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ABSTRACT 

This paper presents infrared spectra taken with the newly commissioned NIRSPEC spectrograph on 
the Keck II Telescope of the High Redshift Radio Galaxy MRC 2025-218 (z=2.63) These observations 
represent the deepest infrared spectra of a radio galaxy to date and have allowed for the detection of 
H/3, [OIII] (4959/5007), [01] (6300), Ha, [Nil] (6548/6583) and [SII] (6716/6713). The Ha emission is 
very broad (FWHM = 9300 km/s) and luminous (2.6x10'*^ ergs/s) and it is very comparable to the line 
widths and strengths of radio loud quasars at the same redshift. This strongly supports AGN unification 
models linking radio galaxies and quasars, although we discuss some of the outstanding differences. The 
[OIII] (5007) line is extremely strong and has extended emission with large relative velocities to the 
nucleus. We also derive that if the extended emission is due to star formation, each knot has a star 
formation rate comparable to a Lyman Break Galaxy at the same redshift. 

Subject headings: galaxies: active — galaxies: structure — galaxies: quasars — galaxies: kinematics 
and dynamics — infrared: galaxies 



1. INTRODUCTION 

Deep radio surveys have proven to be one of the best 
methods for finding high redshift galaxies. Most evidence 
suggests that these powerful radio sources are the precur- 
sors of local giant ellipticals (e.g. Pentericci, et al. 1999). 
Many have irregular and complex morphologies suggestive 
of mergers and they are often surrounded by an overden- 
sity of compact sources; presumably sub-galactic clumps 
(e.g. van Breugel et al. 1998). At both low and high 
redshifts, radio galaxies usually have strong optical emis- 
sion lines, especially OIII at 5007 A. It is strongly debated, 
however, if the emission lines arise by the same mechanism 
as the radio jets. Several authors (e.g. RawHngs & Saun- 
ders 1992, Eales & Rawhngs 1993, and Evans 1998) have 
demonstated a strong correlation between radio luminos- 
ity and [OIII] luminosity, but as Evans showed, there is 
a strong selection effect based on the detection limits as 
a function of distance and this may explain much of the 
correlation. Since the galaxies are often disturbed, star 
formation, large scale shocks and a central AGN are all 
possible sources of the line emission. Active galaxy unifi- 
cation models suggest that radio galaxies are quasars with 
obscured broad line regions (e.g. Antonucci 1993). Eales 
& Rawlings (1993, 1996) and Evans (1998) have been suc- 
cessful at using infrared spectrographs on 4-meter class 
telescopes to measure a few of the brightest lines in small 
samples of radio galaxies in the redshift range 2.2 to 2.6 
and they find line ratios most consistent with Seyfert 2 (ob- 
scured AGN) nuclei. Independent of our current efforts, a 
team has also successfully used the ISAAC instrument on 
the VLT to observe high redshift radio galaxies (HzRGs) 
including MRC 2025-218 (McCarthy, personal communi- 
cation) . 

An additional unexplained phenomenon is that at high 



redshifts (z>0.6) the radio, optical continuum, infrared 
continuum and emission line structures tend to be closely 
aligned (Chambers, Miley, & van Bruegel 1987, and Mc- 
Carthy et al. 1987). This is probably not seen in lower 
redshift targets because the central activity tends to be a 
smaller fraction of the total luminosity than in high red- 
shift sources. Among the proposed explanations are that 
the emission lines arise from shock induced star formation 
(De Young 1989, Rees 1989) or that it is scattered light 
originating from the central nucleus (Fabian 1991). This is 
a crucial question in our understanding of how and when 
most star formation occured in giant elliptical galaxies and 
in clusters in general. 

MRC 2025-218 (z=2.630) has a compact near infrared 
and optical continuum morphology (van Breugel et al. 
1998), but extended Lya emission (5") aligned with its 
radio axis (McCarthy et al. 1992). The extended UV 
emission has significant (8.3±2.3 %) linear polarization 
perpendicular to the UV axis (Cimatti et al. 1996) sug- 
gesting that scattering plays a significant role. McCarthy 
et al. also found three extremely red galaxies (ERO's: R- 
K > 6 mag) within 20" of the radio galaxy. This is a 
large overdensity of such objects and strongly suggests an 
association between the ERO's and the active galaxy. In 
this paper we present infrared spectra taken with a long 
slit oriented close to the radio axis and including one of 
the ERO's. The FRO spectra will be described in a future 
paper. For all calculations we have assumed a cosmology 
with A=0, qo=0.1 and Ho=75 km s"^ Mpc"^ For a red- 
shift of 2.63 this yields a luminosity distance of 2.1x10^ 
Mpc and an angular scale of 7.7 kpc per arcsecond. 

2. OBSERVATIONS AND DATA REDUCTION 

The field of MRC 2025-218 was observed on 4 Jun, 
1999 (UT) with the near infrared spectrograph NIRSPEC 
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(McLean, et al. 1998 and McLean, et al. 2000) on the Keck 
II Telescope during its commissioning. First the field was 
imaged in the K-band with the slitviewing camera which 
is a HgCdTe PICNIC detector (256^ pixels) sensitive from 
1 to 2.5 microns. Figure 1 shows the reduced image of the 
field with a total integration time of 540 seconds and a 
FWHM of 0'.'54. As shown in the figure, the slit (42" long 
and 0"57 wide) was placed on both the radio galaxy and 
the extremely red galaxy dubbed ERO-A by McCarthy et 
al. (1992). This corresponded to a slit position angle of -7 
degrees. 




Aa (arcsec) 

Figure 1. - K band image of the MRC 2025-218 field. 

For spectroscopy, the telescope was repeatedly moved 
roughly 20 arcscconds to center the objects first in the up- 
per portion of the slit then the lower portion. Four 300 
second exposures were taken in both the H-band (~1.6/im) 
and K-band (^2.2fim) yielding an effective integration 
time on MRC 2025-215 of 20 minutes in each band. For 
guiding, NIRSPEC's optical guide camera was used to 
actively track a bright star roughly 2 arcminutes from 
MRC 2025-218. 

Arc lamp and flat lamp spectra were taken at each set- 
ting prior to changing mechanism setups. The 7.6 magni- 
tude AO star PPM 272233 was also observed at the same 
settings in order to remove telluric absorption effects from 
the atmosphere. The calibrator star was reduced first. 
For each band the spectral pair was subtracted and di- 
vided by a reduced flat fleld lamp spectra. Bad pixels were 
then identified and removed by medianing the four near- 
est neighbors. The spectra were spatially rectified using a 
quadratic polynomial at each row, then spectrally rectified 
with a quadratic at each column. The negative spectrum 
of the star was then shifted and subtracted from the posi- 
tive spectrum producing a combined spectrum with resid- 
ual atmospheric lines removed. The stellar spectrum was 
extracted by averaging the central 3 pixels along the 2-d 
spectrum. A synthetic black body spectrum was divided 
into the stellar spectrum and residual hydrogen absorp- 
tion lines from the Brackett series were interpolated over. 
The spectra of the radio galaxy were reduced in a similar 
way except they were divided by the reduced calibration 
star spectrum instead of a black body. For extraction of 
the galaxy spectra, a 6 pixel spatial aperture (1'.'14) was 
used. Spectrophotometry was obtained by determining 



the equivalent widths of the emission lines within a 1'.'5 
aperture in the spectra and comparing this to the broad 
band fluxes of the galaxy in a 1'.'5 circular aperture in the 
slit viewing camera images. 
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Figure 2. - H band spectrum of MRC 2025-218. It is dominated by 
[OIII] at rest wavelength 5007 A. Also present is the other member 
of this doublet ([OIII], 4959 A) and a weak emission line. 



3. RESULTS 

Figure 2 shows the H-band spectrum of MRC 2025-218. 
By far the most dominant line is [OIII] (5007 A) redshifted 
to 1.82 /im. This line is highlighted in figure 3 where the 
complete position velocity map of this line is presented. 
Panel (a) of figure 3 is stretched to highlight the spec- 
trally double nature of the nuclear emission (Av ~ 200 
km sec~^). Panel (b) shows three faint emission knots at 
large angular separations (l"-2") and/or high kinematic 
velocities (~400 km s~^) Although faint, these structures 
repeat in the individual spectra that cover the [OIII] line. 
Two knots appear at essentially km sec~^ relative ve- 
locity, but I'.'S North and 2'.'4 South of the Nucleus. A 
high speed clump appears 1" North of the nucleus and at 
a redshifted relative velocity of 410 km sec~^. This high 
speed clump is also the brightest within our slit with a 
fiux of roughly 1x10""'^^ ergs cm~^. Also detected in 
the H-band spectrum is the other member of the [0111] 
doublet at 4959 A, and H/3. The ratio of [OIII] / H/3 is 
extremely large at 17±7. The H/3 line has a total nuclear 
fiux of only 5 x 10"^'' ergs cm~^ s~^. Table 1 gives all 
detected fluxes and line widths. 

Figure 4 shows the K-band spectrum which is dominated 
by a broad Ha emission line. The spectrum has had a me- 
dian fllter passed over it to improve the appearance of the 
fainter lines. The Ha is well modeled by a pair of Gaus- 
sians having line widths of 9300±900 km/s and 730±100 
km/s. The narrow component is consistent with the Lya 
line width of 700 km/s found by ViUar-Martin et al. 1999. 
After subtracting away the two Ha components, the mid- 
dle graph in figure 4 shows the strong [Nil] (6548/6583 A) 
emission lines as well as weaker features from [01] (6300 
A) and [SII] (6716/6731 A). The line fiuxes and widths are 
also given in table 1. 
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Aslit position (arcsec) 
Figure 3. - Position velocity plots for OIII (5007). Panel (a) is 
stretched to show the double nuclecir peaJc. Panel (b) highlights 
three extended emission regions circled in white. The OIII line is 
highly disturbed with several different kinematic and spatial com- 
ponents including a kinematically split nucleus and a high velocity 
(400 km/s) knot located 2" off nucleus. Nearby OH lines from the 
Earth's atmosphere are labeled. 







TABLE 1 






Emission Line Strengths 






Rest 


Flux 


Line Width 


Line 


A(A) 


(xlO^-*-® ergs s^^ cm~^) 


(km s^"'") 


SII 


6731 


0.4±0.3 


200±100 


SII 


6716 


0.6±0.3 


200±100 


Nil 


6583 


1.3±0.3 


880±100 


Nil 


6548 


1.3±0.3 


880±100 


Ha (narrow) 


6563 


2.7±0.4 


730±100 


Ha (broad) 


6563 


18±2 


9300±900 


OI 


6300 


0.8±0.3 


800±400 


OIII 


5007 


8.4±1.6 


600±200 


OIII 


4959 


2.1±0.4 


600±200 


H/3 


4861 


0.5±0.3 


600±200 



4. DISCUSSION 

4.1. Nuclear Spectrum 

The nuclear spectrum of the HzRG MRC 2025-218 is 
clearly dominated by emission lines from a cxmtral AGN. 
The broad Ha line width is 9300 km/s which is only seen 
in type I AGN (unobscured broad line regions). This line 
width is very close to the mean H/3 line width of 9870±950 
km/s of radio loud quasars in the redshift range 2.0 to 2.5 
by Mcintosh ct al., 1999. The ratio of [0III]/H/3 is 17 
which is also only seen in AGN and ratios of [NII]/Hq and 
[01] /Ha are also consistent with AGN excitation (Oster- 
brock, 1989). 

From the Ha/H/3 narrow line ratio of 5.4 we derive an 
optical extinction Ay=1.4 mag. In this calculation we've 
assumed an intrinsic ratio of Ha/H/3 = 3.1 as seen in local 
AGN (Osterbrock, 1989), and the interstellar extinction 
law of Cardelli et al. (1989). This must be treated as an 
upper limit, however, since radio loud objects may have 
elevated Ha due to coUisional excitation (e.g. Baker et 
al. 1994). If the broad line ratio of Ha/H/3 were simi- 
lar to the narrow line ratio, then broad H/3 should have 
marginally been detected in our H-band spectrum. Wc 
therefore feel safe in the assumption that the extinction to 
the broad line region is similar to the value for the nar- 
row line region (1.4 mag), but not necessarily significantly 
greater. This extinction is also sufficient to explain the 
lack of broad Lya detections in McCarthy et al. (1990) 
and Villar-Martin et al. (1999b). Without extinction our 
broad line emission would predict a Lya broad line flux of 



2.6x10^^'"' ergs s^^ cm^^ in the Villar-Martin slit which 
would have been easily detected but with Ay=1.4 mag 
this is reduced to less than 3xl0~^^ ergs s~^ cm~^ which 
would have been marginably detected at best. 
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Figure 4. - The K band spectrum of MRC 2025-218 is dominated 
by a very wide (9300 km s~^) strong emission line of Ha. The upper 
graph is the reduced spectrum overlayed with the Ha profiles. The 
dashed line is the x-axis for this graph. The middle graph has had 
the continuum and Ha emission lines subtracted to emphasize the 
weaker lines of NIL The bottom graph is the residuals after subtract- 
ing gaussians for each emission line. 

Given the similarities in line width with radio loud 
quasars we now try to determine if the extinction could ex- 
plain the observed differences between MRC 2025-218 and 
radio loud quasars. If we correct the Ha flux for Ai/ = 1.4 
mag (the upper limit to the narrow line extinction), then 
the broad line flux becomes 5.2x10"^^ ergs s~^ cm~^ or a 
broad line Ha luminosity of 2.6x10*^ ergs s~^. We used 
the sample of quasars of Mcintosh et al. (1999) to de- 
rive a mean Ha luminosity of 6.11 xlO**^ ergs s~^ based 
on the mean H/3 equivalent width of their sample, no ex- 
tinction and an intrinsic ratio of 3.1 between Ha and H/3. 
The one sigma dispersion in this value is only 10% in their 
sample. Our extinc;tioii corrected Ha luminosity is then 
weaker than their mean by a factor of 2.6 suggesting the 
central engines are very similar. If we go a step further and 
assume that the intrinsic luminosities are the same, then 
the broad line extinction would need to be Ay=3.5 mag 
instead of Ay=1.4 mag as derived above for the narrow 
line region. 

A remaining difference between MRC 2025-218 and the 
quasars in the Mcintosh sample is the H-band magnitude. 
MRC 2025-218 has a broad band magnitude of H=19.1 
while the mean quasar H-band magnitude is 15.16. After 
correcting for the different redshifts (quasar mean z=2.2) 
then MRC 2025-218 is 3.2 magnitudes fainter than the 
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quasars at a rest wavelength of 4550 A. If the broad band 
flux of MRC 2025-218 is dominated by the AGN then it 
would require Av'=4.2 mag to make it equal to the quasar 
sample. This is surprisingly close to the value of 3.5 mag 
required to match the broad Ha fluxes. The assumption 
that the AGN dominates the broad band flux in a radio 
galaxy, however, is not obvious and may be in conflict with 
the empirically determined K magnitude versus redshift 
relation observed in both low redshift and high redshift 
objects (Bales et al. 1997). MRC 2025-218 is consistent 
with the K vs. Z relation both with and without taking 
the line emission into account. 

Villar-Martin et al. (1999) find that MRC 2025-218 has 
large ratios of [NV]/HeII and [NV]/[CIV] and suggest that 
the most likely explanation is that N is overabundant. 
They held out the possibility, however, that contamina- 
tion from a broad line region was enhancing this line in 
comparison to their other radio galaxies. But they argued 
against this due to the lack of any broad lines including 
[cm]. From our broad Ha detection, however, we clearly 
see that the broad line region is only partially obscured 
and the strong NV emission is probably not indicative of 
high metalicity. This is further corroborated by the rela- 
tively low ratios of [Nil] /Ha (narrow). 

4.2. Spectral Shape and Extended Emission 

The double spectral peak found in [QUI] could be due 
to a high velocity (200 km s~^) cloud of gas or possibly a 
double active nucleus. The unsmoothed Ha narrow line is 
quite noisy but also shows a double profile with a separa- 
tion of 200 km s^i. Due to the noise, however, we are not 
confident in the second Ha peak. If the second peak were 
due to a star forming region it would be unlikely that the 
[OIII] line would be double as well since the OIII/H/3 ratio 
should be much lower for a star burst. 

The off nucleus knots seen in [QUI] are difficult to un- 
derstand. Extended OIII has been observed in other radio 
galaxies aligned to the radio axis (Armus et al. 1998) 
but no line ratios have been determined for this gas. If 
we assume that the emission is from starbursts then our 
brightest knot (1x10^^^ ergs s~^ cm~^) would have an 
[OIII]/Ha ratio less than 1.0. This would make the Ha 
flux greater than lxl0~^® ergs s~^ cm~^ and a luminosity 
more than 5x10"*^ ergs s~^. Assuming the relationship of 
Kennicutt (1983) that the star formation rate is equal to 



L(Ha) divided by 1.12xl0'*i ergs/s we derive a star for- 
mation rate of 45 Mq yr^^. This is comparable to the 
rates seen in Pettini et al. (1998) where they studied 5 
star forming galaxies in the redshift range 2.2 to 3.3. This 
is also close to the estimated star formation rate of the Ly- 
man Break Galaxy MS1512-cB58. As calculated in Teplitz 
et al. (2000) cB58 has a SFR of 620 Mq yr'^ but after re- 
moving a factor of 30 for gravitational lensing this becomes 
21 Mq yr-i. 

5. CONCLUSIONS 

We have obtained the most sensitive infrared spectra 

ever taken of a high redshift radio galaxy. The galaxy 
has very strong emission lines with ratios and line widths 
consistent with an obscured quasar. The narrow line re- 
gion appears to be partially obscured with Ay around 1.4 
mag, but from comparisons with high redshift quasars, we 
estimate that the extinction to the broad line region is be- 
tween 3 and 5 magnitudes. Since other radio galaxies in 
the same redshift range don't show broad emission lines, 
we suggest that MRC 2025-218 is further along in its evo- 
lution towards an unobscured quasar. We cannot rule out 
any of the proposed mechanisms for the production of the 
aligned emission. But based on the [OIII] line strength if 
the majority of the emission is due to star formation, we 
find that the star formation rate is comparable to that of 
Lyman Break Galaxies at similar redshifts. We urge even 
deeper observations of this and other similar radio galaxies 
in order to measure additional extended line emission. 
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